Please note this document is a working draft.

Introduction

Background

Cities have used urban green spaces as a tool to enhance the health of residents, increase biodiversity, and mitigate the effects of climate change. Research has shown that urban green space can alleviate urban-heat-island effects, improve air quality, and encourage outdoor socialization and physical activity.1 However, urban green spaces are more frequently supported in higher-income, predominantly white neighborhoods, exacerbating inequity in health.2 Several U.S. cities, including Denver, have proposed goals to equitably increase green space in their cities.

Often, urban planners and residents use turf and exotic plants to expand green space, which are are not ideal for pollinator health and demand high amounts of water, which may not be sustainable in the mountain west, a region expected to experience more severe drought due to climate change. While a growing body of research suggests exposure to urban green space can improve human health and prevent premature mortality,1,3 little is known about the potential impact of native vegetation on human health in cities in the U.S. mountain west, a region with a growing population vulnerable to the impacts of climate change.

Objective

We aim to estimate the number of premature deaths that would be prevented by the implementation of various policy scenarios adding green space, including native vegetation, to the City of Denver. These policy scenarios were informed by conversations with local stakeholders in advocacy, research, and governance.

Methods

Measurement of existing green space and definition of native-plants

We measured existing green space in the City of Denver by the normalized difference vegetation index (NDVI) measured by the Landsat-8 satellite at a spatial resolution of 30 square meters. To define native-plants greenness, we measured the NDVI of areas in Denver that are mostly native, including the Denver Botanic Gardens and a few other locations on several cloud-free spring and summer days in 2021. Over the course of these days in these areas, the NDVI values tended to range between 0.4 and 0.5 (details here). The NDVI of native vegetation may vary by ecosystem conditions and level of maintenance. To accommodate this uncertainty and assess the sensitivity of results to the NDVI-based definition of native plants, we conducted analyses defining the NDVI of native plants at both 0.4 and 0.5.

Stakeholder engagement to develop policy scenarios

We conducted ten interviews with local stakeholders who are actively engaged in advocacy, policy, and governance related to native plants and green space in Greater Denver. Following these interviews, we developed four policy scenarios informed by the collective set of visions and priorities articulated by the stakeholders. The scenarios range from short-term goals to more ambitious, perhaps less immediately realistic visions. The scenarios fall into five categories.

  1. Homogenous native-plants greening: Under the first set of scenarios, we consider adding native plants homogenously across all block groups, without respect to the specific location. Specifically, we consider two scenarios:

    • Set 30% of the block group’s area to be as green as native plants, corresponding with a goal set by some scientists and advocatess to protect 30% of lands and oceans by 2030[@dinerstein2019];

    • Set 20% of the block group’s area to be as green as native plants, a level that, according to some of our interviewees is sufficient to support pollinator health.

  2. Equity-driven native plants greening: Under the second set of scenarios, we propose adding native plants to the same amount of area as above (30% or 20% of Denver’s area) but rather than applying the intervention homogeneously across all black groups, we prioritize application of the native plants intervention to those areas of greatest need from a health-equity perspective. Both the State of Colorado and the City of Denver have created equity indices to guide their work. The Colorado Department of Public Health & Environment (CDPHE) describes their equity index here4 and makes the index available for download here.5 The City of Denver Department of Public Health & Environment (DDPHE) describes the neighborhood-level equity index here6 and makes the data available for download here.7 We use both definitions to consider a total of four scenarios:

    • Set 30% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the CDPHE definition of health equity.

    • Set 20% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the CDPHE definition of health equity.

    • Set 30% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the DDPHE definition of health equity.

    • Set 20% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the DDPHE definition of health equity.

  3. Create native buffers around riparian areas (rivers, streams, lakes, and ponds) of the following sizes:

    • 200 feet (ideal for ecosystem health but possibly not realistic)

    • 100 feet (compromise)

    • 50 feet (most realistic; less good for ecosystem health)

  4. Initiatives related to green infrastructure and stormwater management. (Please refer to corresponding section for details.)

  5. Replace a portion of parking-lot surface with native plants:

    • 100% of the existing parking area (unrealistic, out of curiosity)
    • 50% of the existing parking area (perhaps realistic)
    • 20% of the existing parking area (most realistic)

Visualization of scenarios

This section elaborates on the scenarios and visualizes them with maps.

Scenario 1: add native plants to all census block groups

Under the first scenario, which we view as the most ambitious, we consider the expected impact on mortality by applying a green-space intervention to all census block groups, without specifying where. This scenario was informed by conversations with colleagues who stated that

We first measured the mean NDVI of each census block group on July 4, 2021. The weighted block-group-level mean is presented here, weighted by the proportion each 30 square-meter pixel covers by the census block group. For example, if half of a pixel overlaps the block group, it receives a weight of 0.5 in the weighted average. We also removed bodies of water before measuring NDVI.

For each definition of native-plants NDVI, We excluded census blocks whose baseline NDVI was above the threshold. This map visualizes the NDVI threshold at 0.4.

And this one defines the native-plants threshold at 0.5.

Scenario 2: add native plants prioritizing census block groups with most health-equity need

Under this set of scenarios, we assumed we had a “budget” of greening of either 30% or 20% of Denver’s area, and we applied that budget to those with most need from a health-equity standpoint based on the CDPHE and DDPHE definitions, described above. As the below maps show, the two health-equity definitions from the State and from the City overlap considerably.

30% of area prioritizing block groups of greatest need according to two health-equity definitions

20% of area prioritizing block groups of greatest need according to two health-equity definitions

Scenario 3: add native-plants buffers to riparian areas

We measured NDVI in a 200-foot buffer, a 100-foot buffer, and a 50-foot buffer around all bodies of water in Denver. We downloaded bodies of water from OpenStreetMap (code here). We defined residential exposure to these riparian areas as those individuals living within a 500-meter buffer, following research on green-space and health that have defined green-space exposure based on residential proximity.3 We estimated the number of people in this buffer by multiplying the population density of the census block group by the intersecting area.

The below map depicts mean NDVI in the portions of census block groups that intersect a 200-foot buffer as well as those pieces that intersect the part of the 500 m buffer that would not be intervened upon, i.e., the part between the 200-foot buffer and the edge of the 500 m buffer.

Scenario 4: green infrastructure

We spoke with representatives at the local Office of Green Infrastructure who described three categories of initiatives–some planned, some aspirational–that could add native or adapted plants to the city.

  • Large stormwater retention projects:
    • Large ponds or basin located on public property that treat collect and treat stormwater after it has been collected in a storm pipe. They are usually vegetated, often with native or adapted plants.
    • Our colleague expected about 75% of the facility’s footprint would have native or adapted plants.
  • Green streets; according to our colleague:
    • Historically, about 2.7 miles of green streets each year, and each mile equates to about 0.15 acres of native or vegetated landscape.
    • Short term goal: increase output to 5.0 miles of green streets; same amount of vegetation per green mile
    • Aspirational goal: output to 5.0 miles; increase vegetated area to 0.75 acres per green mile
  • Proposed stormwater controls on new or re-development:
    • Our colleague stated that as properties develop or redevelop, they may be required to include stormwater runoff control measures to offset negative impacts to flooding and water quality downstream of the site associated with the impervious surfaces added during the development. These stormwater control measures are often green-on-the-ground practices vegetated with native or adapted plants.
    • The requirements may differ by parcel size. Our colleague estimated the following number of sites that may require stormwater control by parcel size:
      • greater than 1 acre: about 100 sites per year
      • 0.5-1.0 acre: about 25 sites per year
      • less than 0.5 acre: 400 sites per year

4.1. Large stormwater retention ponds

We were provided a list of planned projects throughout Denver. Per our conversations, we assume that about 75% of the project’s footprint would consist of native or vegetation. Like for riparian areas, we defined exposure to the projects as those individuals living within a 500 m buffer of the projects. A map of baseline NDVI of the projects themselves and of a 500 m buffer around the projects is below.

4.2. Green streets

Not yet complete.

4.3. Site new or re-development storwmater controls

The approach for estimating the health impact of the redevelopment controls differs from scenarios 2 and 3.1. because we will not know where the re-developments will occur. We thus simulated possible locations. Per our conversations, we anticipate the following number of parcels will be subject to these rules per year:

  • greater than 1 acre: about 100 sites per year
  • 0.5-1.0 acre: about 25 sites per year
  • less than 0.5 acre: 400 sites per year

We gathered data on parcels from Denver’s Open Data Portal (Existing Landuse 2018) and measured their area. A subset of these parcels near Union Station is mapped below.

Then, from all parcels, we sampled 100 sites of size greater than 1 acre, 25 sites of size 0.5-1 acre, and 400 sites less than 0.5 acres. One such sample appears below.

From this point, we followed the same framework as for other scenarios. The baseline NDVI of the sampled parcels and that of their corresponding 500m buffers is visualized below.

Scenario 5: add native plants to parking lots

Finally, we propose replacing a portion of parking-lot surface with native plants. Denver’s area is about 9% parking lot, so even replacing a small amount of parking-lot surface could have a large impact on the total area of urban greening. We obtained spatial data on parking lots](https://www.denvergov.org/opendata/dataset/parking-lots-2016)) from the City of Denver Open Data Catalog.

Total area (mi2), parking lots Total area (mi2), Denver proportion parking lots
13.57 155 0.09

We similarly measured the NDVI on the parking lots as well as the NDVI in the census block groups within a 500-m buffer radius of any parking lot (which is almost the whole city). We present a small subset of these areas near Union Station:

Health-impact assessment methods

We estimated the number of premature deaths averted under each scenario by following a recent meta-analysis that estimated that for every 0.1 unit increase in exposure to NDVI, the relative risk of premature death decreases about 4% (pooled risk ratio of 0.96; 95% confidence interval [CI]: 0.94, 0.97).3 We used this risk ratio to estimate the population attributable fraction8 corresponding to the proposed changes in NDVI from the baseline levels for each scenario to the alternative native-plants level.

Under the first two scenarios, the areal unit is the block group itself, and because we do not specify where, exactly, in the block group the greening intervention would occur, we define NDVI exposure as the average NDVI in the block group and compare that value with the alternative native-plants value. We use 5-year American Community Survey Data (2016-2020) to estimate the population in 5-year age groups in each census block group. For the other scenarios, we draw a 500-m buffer around the intervention area, as shown above. To estimate the population affected, we multiply the population density of the block group in that age group by the area of the block group covered by the intervention area.

To estimate the number of premature deaths prevented in each age group and area, we multiply the population-attributable fraction by the baseline mortality rate by the population size. We restrict analyses to adults aged 30 and above (subject to change) following the age range of many of the cohort studies reviewed3 and exclude northeastern census block groups near the airport (subject to change; perhaps not for all scenarios).

We estimate uncertainty in both the block-group-level population estimates and the relative risk estimate (95% CI: 0.94, 0.97) by resampling from normal distributions and taking the 2.5th and 97.5th percentiles over 1,000 replications. This code has additional detail.

Results

Results of each scenario appear in the table below. The first table describes, for each scenario, the population affected, the cumulative area of the treatment itself (e.g., the area of the retention ponds or parking lots) and the cumulative area of the 500 m residential buffer surrounding the treatment area. Note for Scenarios 1 and 2, the treatment area and the residential-buffer area are the same (the area of the affected census block groups) because we do not say exactly where in the block groups the greening treatment would occur.

Scenario Sub-scenario NDVI, native Area, treatment (mi2) Area, residential buffer (mi2) Pop. affected, est. Pop. affected, 95% LL Pop. affected, 95% UL
all-bg 100-pct 0.4 1,726 1,726 255,199 249,743 260,703
all-bg 100-pct 0.5 2,510 2,510 376,816 370,334 383,152
all-bg 30-pct 0.4 1,726 1,726 255,199 249,568 259,692
all-bg 30-pct 0.5 2,510 2,510 376,816 371,392 383,864
all-bg 20-pct 0.4 1,726 1,726 255,199 248,935 260,491
all-bg 20-pct 0.5 2,510 2,510 376,816 370,134 382,511
all-bg-equity cdphe-def-30-pct 0.4 742 742 94,057 91,070 97,479
all-bg-equity cdphe-def-30-pct 0.5 792 792 104,970 100,942 108,340
all-bg-equity cdphe-def-20-pct 0.4 520 520 56,989 54,221 59,172
all-bg-equity cdphe-def-20-pct 0.5 538 538 61,145 57,971 62,910
all-bg-equity denver-def-30-pct 0.4 697 697 86,153 82,199 88,682
all-bg-equity denver-def-30-pct 0.5 792 792 107,015 101,846 108,952
all-bg-equity denver-def-20-pct 0.4 515 515 54,420 51,558 56,152
all-bg-equity denver-def-20-pct 0.5 538 538 60,046 57,292 62,240
riparian 200-ft 0.4 104 895 113,822 110,362 117,645
riparian 200-ft 0.5 153 1,309 177,520 173,215 181,115
riparian 100-ft 0.4 51 895 113,822 110,826 117,409
riparian 100-ft 0.5 75 1,309 177,520 173,080 181,466
riparian 50-ft 0.4 25 888 113,334 110,285 116,549
riparian 50-ft 0.5 37 1,309 177,520 173,035 181,157
ogi ogi_proj 0.4 14 376 65,527 62,960 67,836
ogi ogi_proj 0.5 17 440 74,882 72,377 77,708
ogi parcel 0.4 2 1,381 213,071 208,841 217,775
ogi parcel 0.5 3 2,030 312,936 307,475 317,817
prkng 100-pct 0.4 22 1,741 256,223 250,784 262,189
prkng 100-pct 0.5 24 2,483 370,940 364,248 377,751
prkng 50-pct 0.4 22 1,741 256,223 250,866 261,382
prkng 50-pct 0.5 24 2,483 370,940 364,428 377,396
prkng 20-pct 0.4 22 1,741 256,223 250,089 261,600
prkng 20-pct 0.5 24 2,483 370,940 364,937 376,581

The second table shows, for each scenario, the estimated number of deaths prevented and the number of deaths prevented per 100,000 people affected.

Scenario Sub-scenario NDVI, native Deaths prevented, est. Deaths prevented, 95% LL Deaths prevented, 95% UL Deaths prevented per 100k, est. Deaths prevented per 100k, 95% LL Deaths prevented per 100k, 95% UL
all-bg 100-pct 0.4 95 84 107 37 33 42
all-bg 100-pct 0.5 234 218 258 62 58 68
all-bg 30-pct 0.4 28 25 31 11 10 12
all-bg 30-pct 0.5 68 63 75 18 17 20
all-bg 20-pct 0.4 19 16 21 7 6 8
all-bg 20-pct 0.5 45 42 49 12 11 13
all-bg-equity cdphe-def-30-pct 0.4 11 9 12 12 10 13
all-bg-equity cdphe-def-30-pct 0.5 24 21 26 22 20 25
all-bg-equity cdphe-def-20-pct 0.4 4 4 5 8 7 9
all-bg-equity cdphe-def-20-pct 0.5 9 8 10 15 13 16
all-bg-equity denver-def-30-pct 0.4 8 7 9 9 8 10
all-bg-equity denver-def-30-pct 0.5 20 18 21 18 17 20
all-bg-equity denver-def-20-pct 0.4 4 3 4 7 6 7
all-bg-equity denver-def-20-pct 0.5 8 7 8 13 11 14
riparian 200-ft 0.4 3 2 3 2 2 3
riparian 200-ft 0.5 8 6 9 4 3 5
riparian 100-ft 0.4 1 0 1 0 0 1
riparian 100-ft 0.5 2 2 3 1 1 2
riparian 50-ft 0.4 0 0 0 0 0 0
riparian 50-ft 0.5 1 1 1 0 0 1
ogi ogi_proj 0.4 1 0 1 1 0 2
ogi ogi_proj 0.5 1 0 2 2 1 3
ogi parcel 0.4 1 1 1 0 0 0
ogi parcel 0.5 1 1 1 0 0 0
prkng 100-pct 0.4 25 22 28 10 9 11
prkng 100-pct 0.5 42 38 46 11 10 12
prkng 50-pct 0.4 12 11 14 5 4 5
prkng 50-pct 0.5 21 19 23 6 5 6
prkng 20-pct 0.4 5 4 5 2 2 2
prkng 20-pct 0.5 8 8 9 2 2 2

Next steps

As of April 27, 2022:

  • Add green streets HIA

  • Add economic-impact attributable to health

References

1.
2.
Wen M, Zhang X, Harris CD, Holt JB, Croft JB. Spatial Disparities in the Distribution of Parks and Green Spaces in the USA. Annals of Behavioral Medicine. 2013;45(S1):18-27. doi:10.1007/s12160-012-9426-x
3.
Rojas-Rueda D, Nieuwenhuijsen MJ, Gascon M, Perez-Leon D, Mudu P. Green spaces and mortality: A systematic review and meta-analysis of cohort studies. The Lancet Planetary Health. 2019;3(11):e469-e477. doi:10.1016/S2542-5196(19)30215-3
4.
CDPHE. Climate Equity Data Viewer. Published online September 13, 2021. https://storymaps.arcgis.com/stories/46bf289f92bc4629a0a1266de4bb7f97
5.
6.
7.
Denver. Denver open data catalog: Equity index 2020 - neighborhood. https://www.denvergov.org/opendata/dataset/city-and-county-of-denver-equity-index-2020-neighborhood
8.
Murray CJ, Ezzati M, Lopez AD, Rodgers A, Vander Hoorn S. Comparative quantification of health risks: Conceptual framework and methodological issues. Population Health Metrics. 2003;1(1):1. doi:10.1186/1478-7954-1-1


Copyright © 2022 Michael D. Garber